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1 Abstract The fluid dynamics of a water jet impinging on a flat surface in confined conditions 

2 were studied using particle image velocimetery. The experiments were meant to replicate 

3 conditions expected in a jet erosion test (JET) designed to assess cohesive sediment ero- 

4 sion parameters in field applications. High-resolution two-dimensional velocity vectors were 

5 measured in a plane passing the jet centerline including free jet, impingement, and wall jet 
e regions within a fixed- wall box. The general flow behavior in the free jet and wall jet regions 
7 is in good agreement with the behavior of impinging jets in an unconfined environment, 
s Results show that the entrainment coefficient, however, is lower than values in unconfined 

9 conditions, lowering lateral spreading rates. The rate of momentum transfer also increases 

10 along the axial direction since the confinement causes secondary flow and recirculation in 

11 the box. Wall shear stress is calculated based on extrapolation of Reynolds shear stress and 

12 turbulent kinetic energy, where the latter procedure provides more consistent results with 

13 expected scour hole shape under an impinging jet. This wall shear stress distribution shows 

14 higher values near jet impingement in comparison to previously reported distributions, espe- 

15 dally that formulated for the JET under unconfined conditions. The maximum value of wall 

16 shear stress is found to be about 2.4 times greater than the commonly accepted value in the 

17 literature, and also occurs at a position closer to the impingement point. The shear stress at the 
is impingement point is also close to its maximum value, which is consistent with the expected 

19 scour hole shape beneath an impinging jet. These findings have important implications for 

20 the use of jet impingement theory to assess sediment erosion, especially in the application 

21 of the JET. 

22 Keywords Impingement jet • Wall shear stress • Confinement • Soil erosion • 

23 Field apparatus 
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1 Introduction 

Impinging jets play an important role in many industrial, mechanical, and environmental engi- 
neering processes such as heating, cooling, drying, rotorcraft brownout, headcut erosion, and 
discharge of pollutants in rivers, lakes, and oceans [1-3]. One of the practical applications 
of an impinging jet is in the assessment of soil erosion and erodibility. Many studies have 
employed an impinging jet to investigate the erodibility of soils and sediments in stream 
channels, dam embankments, and marine and estuarine environments [4-11], although as 
described below, these studies generally are based on results obtained in unconfined condi- 
tions, whereas the apparatus used for testing imposes confinement. 

A submerged jet erosion test (JET) was designed for testing materials in the field [12, 13], 
and is now a developed ASTM standard used commonly around the world [14]. The goal 
of in situ testing devices such as this is to deliver an accurate estimation of soil erodibility 
without disturbance to the material that may occur in moving a soil sample to the laboratory. 
Sample disruption due to transport is a potential problem in other testing devices such as the 
rotating cylinder apparatus [4], flume experiments [15], drill hole apparatus [16], Sedflume 
[17], or the erosion function apparatus [18]. The JET can be applied for either horizontal 
(stream bed) or angled (side channel) surfaces [12,13]. The device uses a constant head 
jet mounted inside an enclosed cylinder that defines the test space. The volume of soil 
removed under the impinging jet (as the scour hole deepens) is determined using a point 
gauge. Based on the scour hole depth, the rate of erosion is determined. Jet impingement 
theory then is applied to determine the shear stress on the soil surface, and the critical 
tractive shear stress and the soil's erodibility coefficient are derived. This device has been 
used in many locations worldwide to quantify the erodibility of soils, providing researchers 
with a tool to evaluate in situ erosion potential [19,20]. Contemporary research, however, 
is still aimed at developing better understanding of the phenomena involved with the JET 



The equations used in conjunction with the JET are based on jet impingement theory, 
which considers an axisymmetric, circular, impinging jet with three distinct regions of flow 
as shown in Fig. 1 [22]. The first is the free jet region (I), where the impingement plane has 
no effect. The second is the impingement region (II), where the jet significantly decelerates 
as it approaches the impingement plane, and then is directed from the axial direction to the 
radial direction. After some radial distance from the impingement point the wall jet region 
(III) develops. 



[21]. 




d 



Fig. 1 Axisymmetric impinging jet and three regions of interest (after [22]) 
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57 Using the JET to assess soil erosion requires an accurate understanding of jet hydrody- 

58 namics. Although many researchers have studied impinging jets, a number of questions arise 

59 in the application of impinging jet theory to assess sediment entrainment. For example, even 
eo when the bed is flat as at the beginning of the test, Phares et al. [23] have shown that the 

61 commonly accepted distribution of applied bed shear stress as determined by Beltaos and 

62 Rajaratnam [22], and assumed by Hanson [12], was derived incorrectly (see below). Also, 

63 the real-time applied shear stress is not constant when the jet is eroding the surface and the 

64 surface is no longer flat (e.g., [24]). Furthermore, a review of impinging jet studies in the 

65 literature reveals that (1) all studies have been done with an air jet except Hanson et al. [25], 
ee (2) no studies have been conducted using the actual JET employed in soil erosion assessment, 
67 and (3) previous work has not considered the possible effects of the confined environment 
es on jet hydrodynamics and the scaling arguments for velocity and impingement shear stress. 

69 This study focuses on this last issue and investigates the flow characteristics of an impinging 

70 jet within a confined space similar to the JET, so that the results can be used to verify the 

71 relations and coefficients employed to quantify sediment erosion indices in the field. The 

72 main objectives of this paper are (1) to characterize the time-mean and turbulent flow of a 

73 circular impinging jet within a confined volume, (2) to quantify the shear stress distribution 

74 on a flat bed surface at impingement and with distance away from this point, and (3) to 

75 critically assess the application of an impinging jet for sediment erosion measurement and 

76 characterization. 



77 2 Experimental methods 

78 Experiments were conducted using a box of size 670 x 670 x 400 mm (4.24/7 x 4.24// x 

79 2.53//), where H = 158 mm is the impingement height). The confinement ratio (box area 
so to nozzle area) was 13,950, and a V-notch weir on the left side of box (viewed from the 
si camera position) provided the outlet for the box. The laboratory setup included a particle 

82 image velocimetery (PIV) system, test box, jet device, injecting syringe, basin, and pump 

83 (Fig. 2). 

84 The conditions for the JET in the field were reproduced in the laboratory using a jet with 

85 a 6.4-mm nozzle diameter do. The test volume here, however, was slightly larger than the 
se JET, and flat vertical surfaces were used (in contrast to the cylindrical volume of the JET) 
87 to facilitate the acquisition of PIV data. Water was pumped from a basin to a constant head 
ss device to produce the jet (see Fig. 2b). A syringe was used to inject a solution of seed particles 
89 into the water just before the nozzle, since the main interest was in the jet flow and not the 
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Table 1 General properties of the jet where Uq is the ideal nozzle velocity (no losses), Uq is the measured 
nozzle velocity, K is the loss coefficient, <2o is m e nozzle discharge, and Mq is the nozzle rate of momentum 
transfer 



Experiment 


h 0 (m) 


d 0 (m) 


^ (m/s) 


Uq (m/s) 


K 


Qo (m 3 /s) 


M 0 (N) 


Run 1 (Low head) 


0.735 


0.0064 


3.80 


3.67 


0.07 


104 x 10 -6 


0.385 


Run 2 (High head) 


1.487 


0.0064 


5.40 


5.27 


0.05 


149 x 10 -6 


0.785 



90 ambient fluid. The experiments were done using two heads ho: a low head (Run 1) of 735 

91 mm, and a high head (Run 2) of 1,487 mm. Using Bernoulli's equation, the ideal (no losses) 

92 nozzle velocities Uq for Run 1 and Run 2 were 3.80 and 5.40 m/s, respectively. These and 

93 other experimental conditions are summarized in Table 1 . 

94 The PIV system was used to measure the two-dimensional velocity field focusing on the 

95 plane passing the jet centerline (Fig. 1), and was measured in five separate zones: (1) the 

96 near-field of the free jet region, (2) the far-field of the free jet region, (3) the impingement 

97 region, (4) the near-field of the wall jet region, and (5) the far-field of the wall jet region. 

98 This procedure allowed greater accuracy and spatial resolution than could be obtained with 

99 a single image of the entire flow. Note that the wall jet results were developed for the right 

100 side of the jet, opposite to the discharge outlet. In addition, small discontinuities sometimes 

101 arose when combining two data sets that overlapped slightly in space (see below). The spatial 

102 resolution of Zones 1 and 2 were 0.9 and 2.1 mm, respectively. The resolution of Zone 1 

103 was chosen to be higher than that of Zone 2 because of high velocity gradients near the 

104 nozzle. A high spatial resolution of 0.6 mm also was chosen for Zones 3, 4, and 5 because of 

105 interest in the wall shear stress. The jet was allowed to run for at least one hydraulic residence 

106 time (test box volume/jet flow) before any measurements were taken, and preliminary tests 

107 confirmed there were no substantial changes in jet characteristics after that time. The camera 

108 was able to capture up to 500 images in memory for each data collection period. Because 

109 of laser warm-up time, variation in particle injection, and other experimental issues, about 
no 350-450 images were considered valid for each run. The trigger rate of the camera was 
m set to 80 Hz, and 13 data collection periods for each zone were captured for averaging 

112 purposes. Out of those periods, 10 or 11 data collection periods in Zones 1 and 2, and 5 

113 or 6 data collection periods in Zones 3, 4, and 5 typically were identified as valid. The 

114 smaller number of valid images for measurements close to the wall is thought to be due to 

115 interferences related to the presence of the wall. The total lengths of time of captured flow 

116 used in the analyses for the free jet and wall jet regions were thus 60 and 30 s, respectively. 

117 Exposure time for each laser frame and the time between pulses were set at 120 and 100 jis, 
1 1 s respectively. Dynamic Studio Version 13.4 software was used to analyze images . The software 

119 was able to process images and build 2D- vector maps in the domain. Spatial calibration was 

120 accomplished using a target of 270 x 190 mm with white dots and black background, and 

121 an interrogation area of 16 x 16 pixels was chosen for all zones except Zones 1 and 2 

122 of Run 2 which have an interrogation area of 32 x 32 pixels. Based on Willert [26], the 

123 uncertainty of PIV displacement measurements for in-plane components is on the order 

124 of 0.06 pixel, which is equivalent to 0.28 mm/s (0.56 mm/s) for Zone 1, 0.63 mm/s (1.26 

125 mm/s) for Zone 2, and 0.18 mm/s (0.18 mm/s) for Zones 3, 4, and 5 of Run 1 (Run 2). 

126 Processed data were exported as text files and analyzed using user-defined codes in MATLAB 

127 R2012b. 
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128 3 Results 

129 To characterize jet behavior, the jet properties are first compared to those reported in the 

130 literature for the free jet, impingement, and wall jet regions. Note that do is the only scaling 

131 length for the free jet region, while for the impingement and wall jet regions do and H are 

132 alternative scaling lengths, where H/do = 24.1 for the present tests. For the following plots 

133 that have alternative distance axes, the bottom axis is the distance normalized using the 

134 default scaling length for the plots, and the top axis is an additional scaled distance to better 

135 understand the relation between free and wall jet regions. 

136 3.1 Free j et (Region I) 

137 The free jet definition applies to neutrally buoyant turbulent jets diffusing in an unbounded 

138 environment (e.g. [27]). Although the present case is an impinging jet in a bounded environ- 

139 ment, results were compared to the behavior of free jets to identify any differences between 
mo unconfined free jets and the present jet. A submerged free jet has two distinct sub-regions 

141 [28] . The first sub-region is the zone of flow establishment (ZFE), also known as the potential 

142 core, and extends from the jet origin (/ = 0, where / is axial distance from the nozzle) to a 

143 point along the jet axis where the mixing region on all sides of the jet has penetrated to the 

144 centerline of the jet (J = J p , where J p is the length of the potential core, see Fig. 1). Hanson 

145 et al. [25] found J p = 6. 3 Jo- In this region, the time-averaged centerline velocity u m , which 

146 also is the local maximum axial velocity in free jets, is constant and equal to the nozzle exit 

147 velocity Uo- The other sub-region is the zone of established flow (ZEF), which for impinging 
us jets extends to a point where the boundary begins to affect the velocity magnitude and direc- 

149 tion. In this region, u m decreases with increasing /, and the cross-sectional velocity profiles 

150 have a Gaussian distribution. Typically, the ZEF extends to J/H = 0.86 [22]. In order to 

151 obtain high-resolution data, free jet measurements were divided into near-field and far-field 

152 regions and separate measurements were performed for each. Any apparent discontinuities 

153 in plots at J/ do = 8.7 are due to joining those datasets. 

154 3.1.1 General observations 

155 The nozzle velocities Uo are given in Table 1. These measured velocities are just 3.4 and 

156 2.5 %, respectively, lower than the theoretical computed velocities assuming no losses. Thus, 

157 the head loss coefficients K for Run 1 and Run 2 are 0.07 and 0.05, respectively. These 

158 values lead to a discharge coefficient C for the nozzle (Uo = Cjd^^Jlgho) of about 0.97 

159 for both runs, which confirms the range of 0.95 to 1 .00 reported by Al-Madhhachi et al. [29] 

160 for the JET. Measured Uo is used to scale other parameters in the free jet and wall jet regions. 

161 Comparing Run 1 and Run 2, both free-jets evolved in almost the same way with increasing 

162 distance from the nozzle. 

163 3.1.2 Length of the potential core 

164 The length of the ZFE was used by Hanson et al. [20] to relate the maximum wall shear stress 

165 to the maximum stress at the nozzle, so its value is important in how data are presented, as 

166 discussed further below. This length was found to be J p = 2. 3 Jo for both runs (compare 

167 with Hanson's [24] value of 6.3do), and even at J /do = 5 the variation of u m from Uo is 

168 only about 10 %. 
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169 3.1.3 Centerline velocity decay 



170 In the ZEF, the mean centerline velocity decays as [30-32] 




-l 



171 



(1) 



172 where Q is the diffusion coefficient and Jo is a virtual origin for the velocity profile of the 

173 ZEF. Both Cd and To depend on the exit conditions [33]. The reported range for Cj is 5.8 to 

174 7.4, with a commonly accepted average of 6.3 [22]. Hussein et al. [31] proposed Jo = Ado. 

175 Many studies, however, have suggested that the virtual origin may be taken as zero because 

176 it is much smaller than J [28]. For a highly confined free jet in a pipe, Liu et al. [34] found 

177 a slightly modified form of Eq. (1) was applicable, 



179 where D is the diameter of the confining pipe. In this study, D is defined as the equivalent 

180 diameter of the confining box, so that do/D is 0.008. 



182 agreement with Eq. (1) for distances greater than J /do = 4.6 (Fig. 3; note that the vertical 

183 axis is inverted from Eq. 2, to be consistent with previous studies). The values of Q and Jo 

184 in Eq. (1) for the present study, however, are not the same as suggested for unconfined jets. 

185 The free jet is fully developed before entering the impingement region and is affected by 

186 the impingement plane at J = 21.2do = 0.860// for both runs. The fitted decay line was 

187 calculated using the centerline velocity in the ZEF for J /do = 4.6-21.2 for both runs. The 

188 virtual origins of the velocity profile were found to be Jo /do = —2.8 and Jo /do = —2.9 for 

189 Run 1 and Run 2, respectively - although Jo is assumed as 0 below, these values are reported 

190 here to make the data reproducible for future studies. The diffusion coefficient Q was found 

191 to be 7.9 and 7.6 for Run 1 and Run 2, respectively, somewhat greater than in the range of 

192 5.8-7.4 noted above. Figure 3 shows that for the conditions of this apparatus, the centerline 

193 velocity starts to decay in a shorter distance from the nozzle than for ideal free jets. These 

194 differences are believed to be a result of confinement on the jet as discussed below. 

195 Assuming Jo = 0 in Eq. (1) and u m = Uo at / = J p , Hanson et al. [20] wrote 

196 J p = Cddo (3) 

197 As can be seen in Fig. 3, although Eq. (1) predicts the centerline velocity in the ZEF, it is not 

198 valid close to the potential core, and there is a significant difference between measured and 

199 predicted velocities dXJ = J p . Actually, there is a transition zone (TZ) between the ZFE and 

200 ZEF (J /do = 2.3-4.6 in the present case; see Fig. 3) where characteristics of neither zone 

201 apply. The centerline velocity decays by about 9 % in the TZ. By ignoring this transition 

202 length, Eq. (3) was suggested by Hanson et al. [20] to determine J p , while it is more accurate 

203 to say this equation gives the total length of the ZFE and the TZ. Note that although J p was 

204 found to be 2.3 Jo in the present case, it could be a function of parameters such as confinement 

205 ratio and nozzle properties. 

206 3.1.4 Lateral spreading rate 

207 As the distance from the nozzle increases, more ambient fluid is entrained across the jet 

208 boundary, and the jet width increases. The jet nominal boundaries are usually specified by 



178 




(2) 



181 



The present centerline velocity data do not agree with (Eq. 2), although there is good 
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0.8 



+ Run 1 
Run 2 

---Eq. 1 fitted to Run 1 data 
Eq. 1 fitted to Run 2 data 
— Eq. 1 - C d =6.3 and J 0 =4d 0 
— Eq. 2 - C d =6.3 




Fig. 3 Centerline velocity decay as a function of distance from the nozzle for Regions I and II {vertical dashed 
line indicates the boundary between Regions I and II, J/H = 0.86). The zone of flow establishment (ZFE) 
and the transition zone (TZ) are shown on the figure 



209 the location where the time-mean axial velocity u is an arbitrary fraction of the centerline 

210 velocity (l/e, or more usually 1/2, [35]). Lateral spreading in unconfined free turbulent jets 

211 has been shown to be approximately linear and can be expressed as [36] 

r _m = Cu (4 ) 

do U \ do J 

213 where r\/2 is the jet half- width, defined by the location where the time-mean axial velocity 

214 is 1/2 of the time-mean centerline velocity, C u = dri/2/dJ is the spreading rate, and is 

215 the virtual origin for the lateral spreading rate. Values for C u have found to be 0.107 [37] or 

216 0.105 [38]. 

217 The free-jet half- width is shown in Fig. 4 for the two experimental runs. Although the 

218 lateral spreading of neither run is purely linear, a best fit line through the data in the ZEF 

219 (J I do = 4.6 to 21.2) gives values of 0.066 and 0.076 for C u , and values of —2.8 and —1.6 

220 for J^/do in Run 1 and Run 2, respectively. The average spreading rate for the two runs is 

221 about 33 % less than the measurements by Papanicolaou and List [37] and Hongwei [38], 

222 indicating less entrainment as the jet penetrates into the medium, as discussed further below. 

223 3.1.5 Axial velocity distribution 

224 A self- similarity solution to the equation of motion, which also has been shown experimen- 

225 tally, gives the radial distribution of axial velocity for axisymmetric free jets in unconfined 

226 conditions as [22] 

^ = e - a693 "» (5) 
u m 

228 where rj u = r/r\/2 and r is the radial distance from the jet centerline. Figure 5 shows that 

229 the present results are in good agreement with the self- similarity profile not only in the ZEF 

230 but also in the impingement region very close to the bed, although there is some additional 

231 spreading for the profile closest to the bed (J/H = 0.95). 

Springer 



Journal: 10652-EFMC Article No.: 9354 □ TYPESET □ DISK □ LE □ CP Disp.:2014/4/28 Pages: 24 Layout: Small 



Environ Fluid Mech 



J/H 




10 15 
J/d 0 

Fig. 4 Free jet half-width for Run 1 and 2 compared to that in unconfined conditions (vertical dashed line 
indicates the boundary between Regions I and II, J/H = 0.86). The zone of flow establishment (ZFE) and 
the transition zone (TZ) are shown on the figure 
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° J = 0.20H 

v j = oaoih 

o j = o.60# 

« J = 0.80F 
a J= 0.95H. 



0 



Fig. 5 Free jet self- similarity conditions (Run 1: open symbols; Run 2: filled symbols) 



232 3.2 Impingement region (Region II) 

233 The impingement region of a circular jet has been studied extensively [39-43] . In this region, 

234 the impingement plane affects the jet flow, and the flow becomes parallel to the surface at the 

235 end of this region. For large impingement heights H > 8.3 Jo [44], the flow is fully developed 

236 before entering the impingement region, and the critical flow parameters are Uo and do. The 

237 characteristic length scale for the impingement region length, and the pressure and shear 

238 stress distributions on the wall in this region is H [22,45]. Centerline velocity data (Fig. 3) 

239 indicate that the wall affects the flow at J/H = 0.86, which is taken here as the beginning 

240 of the impingement region. 

241 Reynolds stresses, including shear stress u'v' and normal stresses u a and i/ 2 , in the axial 

242 and radial directions, respectively, were not affected by the wall as much as the mean velocity 

243 at the beginning of the impingement region. To assess the impact of the impingement wall 
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A - Run 2 
* v' 2 - Run 2 



ttV - Run^ . 



o 




).5 



0.6 



0.7 



0.8 



0.9 



1 



Fig. 6 Local maximum turbulence stresses scaled by their respective values at J/H = 0.86 

244 on the Reynolds stresses, local maximum stresses u' 2 m , v' 2 m , and u'v' m were normalized 

245 by their respective values at J/H = 0.86. All three stress scales decrease at relatively low 

246 rates as the jet diffuses from the nozzle to the impingement plane, but u' 2 and u'v' decline 

247 rapidly in the vicinity of the wall (Fig. 6). As streamlines turn from the vertical to horizontal 

248 direction in the impingement region, v' 2 increases close to the boundary. The wall affects 

249 u' 2 , v' 2 , and u'v' at J/H = 0.96, 0.91 and 0.92, respectively, causing decreasing u' 2 and u'v' 

250 and increasing v' 2 for larger J/H. For unconfined conditions the stresses are affected by the 

251 wall at J/H = 0.96 [41]. 

252 3.3 Wall jet region (Region III) 

253 The impingement region ends and the wall jet region starts where the excess stress on the 

254 wall approaches zero and the static pressure becomes equal to the ambient pressure. This 

255 transition occurs dXr/H = 0.22 in an unconfined environment [22]. Wall jet flow depends 

256 on H even at large radial distances from the impingement point for large impinging heights 

257 (H > 8. 3 Jo) [46]. Similar to the above, any discontinuities in plots at r/H = 0.4 are due to 

258 joining two datasets collected for this region in post-processing. 

259 3.3.1 Wall jet velocity distribution 

260 It has been shown that the velocity distribution in the wall jet region can be described by a 

261 self- similarity function [35,47], 



263 where v m is the (time averaged) local maximum of radial velocity, z is the vertical distance 

264 from the impingement wall, and the length scale 8 is the jet half- width, which is taken as 

265 the value of z where the radial velocity is 0.5i3 m (note that there are two locations where 

266 v = 0.5t> m ; here the farther one from the wall is considered). Measured radial velocity 

267 distributions follow this self- similarity profile in the wall jet region, verifying the overall 

268 agreement to wall jet behavior as indicated by Eq. (6) (Fig. 7a). The self- similarity profile is 

269 valid for r/H > 0.15. 



262 




(6) 
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270 3.3.2 Velocity and length scales 

271 The length scale 8 [46] and the velocity scale v m [35] are defined for an unbounded wall jet 

272 as 



8 /r\0.9 
— = 0.09f 
H 

V m 1.03 



/ r \ w 

-= 0.098 (-) (7) 

(8) 



Uo r/do 

275 The present results show that the minimum value for 8 (equal to 0.61 do) occurs at r/H = 0.22, 

276 and it increases for greater r (Fig. 7b). Although it follows the nearly linear trend of Eq. 7 

277 for r/H > 0.31, the values are smaller and shifted to the right (Fig. 7b), showing that the 

278 wall jet in confined conditions is not as wide as in unconfined conditions. To account for this 

279 shift, Eq. (7) can be modified by adding a virtual origin as 

8 / r \0.9 

-= 0.098 (-) -no (9) 

281 where ro is a non-dimensional virtual origin and here is found to be 0.007. In a confined 

282 volume, the entrainment is limited by the boundaries and the existence of secondary flow. 

283 This lower entrainment rate results in a lower lateral spreading and a thinner jet. The wall jet 

284 velocity scale follows the overall behavior described by Eq. (8). Its maximum value occurs at 

285 r/H = 0.14, where the wall jet is not yet developed (Fig. 7c). As previously noted, Beltaos 

286 and Rajaratnam [22] proposed that the wall jet region starts dXr/H = 0.22. Although the 

287 minimum jet half- width occurred dXr/H = 0.22, the wall jet is not fully developed in the 

288 present study until r/H = 0.27, where the velocity normal to the wall becomes zero. 

289 3.3.3 Reynolds shear stress behavior near wall 

290 The maximum magnitudes of Reynolds shear stress occur close to the wall (Fig. 8a). The 

291 location of the local maximum occurs at 0.24do above the bed near the impingement point, 
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Fig. 8 Reynolds shear stress characteristics for Run 1: a scaled Reynolds shear stress profiles (length scale 
^mV * s ^ e J et half- width based on u'v' profiles, taken as the value of z where u'v' is 0.5z/ v' m ), b location of 
local maximum (distance from the wall), c local maximum value of Reynolds stress, and d measured Reynolds 
stress at the closest measurement point to the wall (0.4 mm above the wall) 



292 and then moves farther from the wall with greater radial distance (Fig. 8b). Even far from 

293 the impingement point, the local maxima occur at distances less than 2do from the wall. The 

294 stresses are negative by convention. The magnitude of these local maxima increases with r, 

295 and reaches a peak at r/H = 0.26 (Fig. 8c). The measured Reynolds stresses at the closest 

296 position to the wall (0.4 mm above the wall) are shown in Fig. 8d. Near the impingement 

297 point up to a radial distance of 0.40//, the measured values vary between —1.5 and 3.5 Pa, 

298 assuming water density is 1000 kg/m 3 . This large range of values casts doubt on the utility 

299 of using the near-bed Reynolds stress measurements to find the total wall shear stress. 

300 4 Discussion 

301 Results demonstrate that the free jet and wall jet flow of the JET follow the overall behaviors 

302 observed in previous studies for air jets in unconfined environments, but several important 

303 differences are found. Transport of momentum and variation of mass flux along the jet are 

304 considered in this section, followed by a discussion of the distribution of the wall shear stress 

305 and its application to soil erosion characterization. 

306 4.1 Momentum and entrainment 

307 The nozzle discharge Qo and rate of momentum transfer Mq of the jet can be found based on 

308 the nozzle velocity Uo and diameter do (see Table 1). As the jet moves into the box, the flow 
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309 discharge increases because of entrainment of ambient fluid. The entrainment rate is related 

310 to the centerline velocity by an entrainment coefficient a, 



312 where b' is the jet full- width [48], defined in this study by assuming the jet extends to 

313 the location where the axial velocity is 1.5 % of centerline velocity. This percentage is 

314 chosen to include as much of the actual jet momentum as possible. Many studies have 

315 provided measured or calculated values of the entrainment coefficient. For a free jet entering 

316 an unconfined volume, proposed values vary between 0.065 and 0.080 [31,49-51]. 

317 The JET in field applications confines the jet in a cylinder, so the jet enters a finite volume 

318 and excess fluid discharges at the top of the cylinder [13]. This setup introduces two issues 

319 not present in unconfined tests: (1) recirculation of fluid because of the enclosure, and (2) 

320 flows caused by fluid exiting the experimental volume. The one main difference between 

321 the present test setup and the JET is that the experimental environment employed here is 

322 rectangular to avoid problems using PIV on a curved surface, whereas the JET is cylindrical. 

323 The box also encloses a larger volume, so any effect of confinement seen here is likely to be 

324 stronger in the JET. 

325 To study the effect of confinement and the outlet condition, the right and left sides of 

326 the jet for Run 1 are considered separately. The velocity vector map reveals the existence 

327 of strong secondary flows in the box (Fig. 9), which has been observed and discussed in 

328 previous studies on confined jets [34,52-54]. This plot also shows there is an effect of the 

329 outlet location as indicated by the asymmetric ambient flow circulation. However, details 

330 of this effect are beyond the scope of the present study, and in any case the ambient flow 

331 is not fully resolved here. Any observed differences between left and right sides of the jet 

332 can be attributed to the effects of the outlet location. The confining walls of the box cause 

333 the wall jet to turn upward, which drives circulation of flow from the wall jet region back 

334 to the free jet region. This turning begins at r/H = 0.27 on the right side of the wall jet, 

335 and r/H = 0.18 on the left side of the wall jet (Fig. 9). The more intensive effect of the 

336 confinement on the left side is associated with outflow conditions. Using PIV images and 

337 velocity vector maps, the jet axis, defined as the loci of local maximum of axial velocity, was 

338 verified to be orthogonal to the wall boundary to ensure that this asymmetric behavior was 

339 not due to any unintended physical deviation of the jet axis from vertical. The jet width on 

340 either side of the jet axis increases with the same slope at first, but gradually the growth rate 

341 for the left-hand side decreases and deviates from the linear trend (Fig. 10). The right side, 

342 however, extends linearly to the point where the jet enters the impingement region, and the 

343 deviation can be seen after J/H = 0.87 (Fig. 10). 

344 As the jet enters the box, entrainment linearly increases the jet flow rate from the initial 

345 value at the nozzle up to a maximum near the boundary of the free-jet and impingement 

346 region at J/H = 0.89 and 0.82 for the right and left side, respectively. The initial flow rate 

347 for both sides is the same and increases at the same rate up to about / = 5do, where the 

348 effect of confinement and other conditions such as the position of the outlet differentiate the 

349 right and left side of the free jet (Fig. 1 1). It can be seen that the flow rate on the left-side is 

350 less than on the right-side. Figure 12 shows the variation of a along the jet calculated using 

351 two methods: (1) using Eq. 10 with Q calculated by the integration of axial velocities over 

352 the respective side area, which is determined by the jet full width, and (2) using a = u e /u m 

353 [35] with measured entrainment velocity u e , taken as the horizontal velocity component at 

354 the edge of the jet where r = b' . The entrainment coefficient of the jet according to the 

355 second method is more than 85 % less than the above-mentioned values for non-confined 




(inb'u 



(10) 



311 



Springer 



3l Journal: 10652-EFMC Article No.: 9354 □ TYPESET □ DISK □ LE □ CP Disp.:2014/4/28 Pages: 24 Layout: Small 



Environ Fluid Mech 




Fig. 9 Two-dimensional velocity field for Run 1 by joining measurements in Zones 1 and 2. Nozzle is located 
at r = 0 and / = 0. Note that the test box edges are at r = 2.12// although the vectors are affected even 
before r = H 
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Fig. 10 Free jet full- width for the left and right sides for Run 1 , calculated assuming the jet boundary is where 
the axial velocity is 1.5 % of centerline velocity 



356 conditions (Fig. 12). Confinement limits the available volume of fluid to entrain and decreases 

357 the spreading rate of the jet. These limited entrainment rates and consequently low jet width 

358 also were observed in previous studies of confined jets [55-57]. 

359 The rate of momentum transfer M on each side of the jet can be determined by integrating 

360 momentum flux over each (left and right sides) respective cross-section area, 

r=b' 

361 M= J mp(u 2 + u a ^)dr (11) 

r=0 

362 where the total momentum transfer rate is the sum of Eq. (11) evaluated for the left and 

363 right sides of the jet. In the case of momentum-conserving jets, M at any distance from the 

364 nozzle must equal the value at the nozzle in the absence of imposed pressure gradients, and 
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Fig. 11 Free jet flow rate for the left and right sides for Run 1, calculated using jet full- width 
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Fig. 12 Free jet entrainment coefficient (a) for the left and right sides for Run 1 (Method- 1 : based on Eq. 10 
using calculated flow rate; Method-2: based on a = u e /u m using measured entrainment velocity at the jet 
boundary r = b r ) 



365 the entrainment that occurs in the radial direction does not change the axial momentum [58]. 

366 Figure 1 3 shows that the rate of axial momentum transfer varies along the jet axis. This effect 

367 is seen on both sides of the jet, but with different magnitudes. Initially, M is the same for both 

368 sides and close to 0.5Mo = 0A9N. Also, M on both sides increases from the nozzle up to 

369 J = S.ldo = 0.33 H, and then decreases slightly up to J = 0.85//, where it is substantially 

370 affected by the impingement plane. The left- side M, however, is less than the right-side 

371 value except when J < 3 Jo- At first glance it appears unusual that M should increase with 

372 increasing distance from the nozzle, but this increase can be explained by considering the 

373 recirculation in the box. In other words, the entrained flow does not enter the jet with zero 

374 axial velocity; rather, it adds to the axial momentum of the free jet. It is obvious that high 

375 pressure gradients close to the impingement plane are responsible for the decline of M in that 
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Fig. 13 Rate of momentum transfer for the left and right sides of the free jet for Run 1 



376 region. As there are many unknowns such as pressure gradients, gain of momentum across 

377 the jet boundary, and its loss, performing a complete momentum balance through the jet is not 

378 possible. It is generally accepted, however, that the jet momentum in the axial direction is not 

379 conserved through the whole length of the jet [59-61], especially in the confined conditions 

380 due to recirculation effects [62]. 

381 These observations of relatively low jet width and high momentum transfer rates repre- 

382 sent significant effects of confinement on the jet behavior. These results demonstrate that the 

383 recirculation and return flow in a confined environment increase the momentum and decrease 

384 the free jet width. These are consequences of entrainment provided by the recirculated flow 

385 instead of being in an infinite environment in which recirculation does not exist [63]. Exper- 

386 imental measurements of Karimipanah [48] showed that when the ratio of the room area to 

387 the nozzle area was larger than 10,400, the jet momentum was close to that of the unconfined 

388 jet. In the present study, this ratio is somewhat larger, around 13,950, but the confinement 

389 still affected the jet momentum. This could be due to the differences in experimental setups 

390 (Karimipanah considered a compressed air jet and two ceiling outlets). 

391 4.2 Shear stress distribution on the impingement plane 



392 
393 



395 
396 



398 
399 



Based on a solution to the equation of motion in the axial direction, Beltaos and Rajaratnam 
[22] suggested the shear stress distribution at the impingement plane is 



r 

tin 



= 0.18 



/l_ e -H4A 2 \ 
(— -)- 9A3ke 



(12) 



is the 



where X = r/H is the dimensionless distance from the impingement point, and r„ 
maximum shear stress given by 

r m =0.16 H 0 = (13) 
(H/d 0 ) 2 

and r m occurs at A, = 0.14. Examining a water jet impinging on a smooth planar surface, 
Hanson et al. [25] found the wall shear stress distribution to be similar to that of an air jet, with 
shear stress also given by Eqs. (12) and (13). Hanson and Cook [13] assumed the maximum 
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401 
402 



shear stress on the bed causes the maximum scour depth beneath the JET, and defined the 
maximum shear stress, 



c f c d 



(Js/doY 



(14) 



where J s is the distance from the nozzle to the soil surface at the jet centerline, which is 
initially J s = H and increases to J s = H e at equilibrium where r m is equal to critical shear 
stress r c , Cf is the coefficient of friction equal to 0.00416, and is 6.3, resulting in Cf Cj 
= 0.16 (see Eq. 13). Values of Cf and Cd were determined based on direct measurements of 
shear stress on a smooth surface beneath an impinging jet using a hot-film probe by Hanson 
et al. [25]. In this study, C s = Cf Cj is defined as the maximum shear stress coefficient for 
easier comparison of maximum wall shear stress. 

The accuracy of any sediment erodibility assessment using the JET is highly dependent on 
the accuracy of estimation of applied shear stress by the jet on the sediment surface (Eq. 14). 
This equation has been developed for a jet impinging on a smooth, flat surface. The scour hole 
formation in the field results in a curved, rough surface. In addition, the assumption of the 
maximum scour hole caused by the maximum shear stress is not consistent with the off-axis 
location of maximum shear stress in Eq. (12). In addition, although Eq. (12) appears to be 
in good agreement with experimental data from an air jet impinging on a solid surface (e.g., 
[22,39]), Phares et al. [23] pointed out an inconsistency in its derivation. The equation of 
motion was solved in the axial direction for the impingement region, but the wall shear stress, 
which is in the radial direction, was mistakenly defined based on that solution. It appears 
serendipitous then that the maximum shear stress relation (Eq. 13) as an empirical relation 
based on the stated experimental results holds within 15 % of the observed data. Those data 
were measured using a Preston tube with an outside diameter of 1.2 mm calibrated in fully 
developed turbulent shear flows, and Phares et al. [23] explained that the severe pressure 
gradients in the impingement region laminarize the boundary layer and so the Preston tube 
cannot be used to measure shear stress within the impingement region. 

Solving the boundary layer equation in the radial direction along the surface, Phares et 
al. [23] derived a model to find the wall shear stress for the laminar boundary layer of an 
axisymmetric fully-developed jet as 



pU. 



T Rp 1/2 
^Ke 0 



(I) =g2 ^) 



(15) 



431 where Reo = Uodo/v is the Reynolds number based on jet discharge values, v is kinematic 

432 viscosity, and g2 is a universal function. Based on this model, the magnitude r m and location 

433 r m of the shear stress peak are 



= 44.6p^Re-^ (|)" 2 , | = o.09 



(16) 



435 For the turbulent boundary layer, Poreh et al. [46] empirically developed the wall shear stress 

436 distribution as 

r _ 1/9 / H\ 2 i/w r \-2.3 



438 Figure 14 shows that the measured turbulence intensities 



near the wall are 



439 relatively high with a peak at the impingement point. Based on the law of the wall [64], there 

440 is a very thin viscous sublayer very close to the wall where the velocity distribution is linear 
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Fig. 14 Turbulence intensity (a) at jet centerline along the free jet axis, and (b) in the vicinity of the impinge- 
ment wall 



441 and the viscous shear stress can be computed based on the slope of this line. Results of the 

442 wall jet velocity distribution (Fig. 7a) show that this layer is very thin, such that the local 

443 maximum radial velocities occur at the second closest measurement position to the wall (1 

444 mm above the wall) and the velocities at the closest position to the wall (0.4 mm above the 

445 wall) are, on average, 53 % of the local maximum value until r/H = 0.27. It should be noted 

446 that measurements within this layer are difficult to obtain. In this study, the wall shear stress 

447 is determined using two methods. 

448 First, Reynolds shear stress (pu f v f ) was computed above the wall and linearly extrapolated 

449 to the wall. As discussed above, measured near- wall values of pu'v' did not show a reasonable 

450 trend and exhibited oscillations between negative and positive values (Fig. 8d). To avoid any 

451 errors in data very close to the bed, the extrapolation was done using two reliable points. The 

452 first point is the location of local maximum value of absolute Reynolds stress. Note that the 

453 Reynolds stress at this point is negative by convention. This local maximum occurs 1.5 mm 

454 above the bed in the vicinity of the impingement point, and then this distance increases almost 

455 linearly, reaching 1 1 mm at r/H = 0.82 (Fig. 8a). The second point is the location of local 

456 maximum radial velocity, where by definition the total shear stress is zero. Thus, the line 
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Fig. 15 Dimensionless wall shear stress based on extrapolation of Reynolds stress and TKE (dashed lineio 
the left is the impingement region and to the right is the wall jet region) 

457 that passes through these two points is projected to the bed to find the total wall shear stress. 

458 For the second method, Haehnel and Dade [2] showed that for an axisymmetric impinging 

459 jet, u'v' = 0.2TKE where TKE = 0.5 (u' 2 + v' 2 + w' 2 ^ is turbulent kinetic energy per unit 

460 mass. Here, two-dimensional velocities were measured on a plane passing through the jet 

461 axis, so shear stress at the bed was determined using TKE = 0.5 \ u' 2 + i/ 2 ), and it was 



462 assumed that the angular component of velocity fluctuations (v w' 2 ) is negligible. Thus, 

463 TKE values were extrapolated linearly to the bed to estimate wall shear stress. 

464 Both of these methods give almost the same distributions of wall shear stress for r/H > 0.5 

465 (Fig. 15). These distributions, however, differ markedly for r/H < 0.5, although the maxima 

466 are almost the same. Results of the first method depend highly on the accuracy of the location 

467 of the local maximum velocity, but there is some uncertainty in its accuracy because of the 

468 resolution of measurements close to the impingement plane. As noted, the local maximum 

469 velocity occurs less than 1 mm from the wall in the region from the impingement point to 

470 r/H = 0.24, while the data resolution is 0.6 mm in this region. The second method, however, 

471 extrapolates the TKE using four data points near the wall other than the one closest to the 

472 wall to determine the corresponding values on the wall. The accuracy of this method is higher 

473 than for the first method as more points are used for the extrapolation and the accuracy of 

474 extrapolation is not highly dependent on data resolution. 

475 While most of the proposed distributions under an impinging jet show zero wall shear 

476 stress at the impingement point, the literature indicates that the erosion rate at the jet centerline 

477 is, in fact, not zero. Zero shear stress on a surface implies no sediment erosion. Thus, the 

478 shear stress distribution of Beltaos and Rajartnam ([22]; Eq. 12) makes little physical sense 

479 close to the jet centerline, where on a soil bed the maximum scour hole depth forms (e.g., 

480 [13,65]). The calculated distributions of shear stress for both Runs 1 and 2 show a non-zero 

481 value at r/H = 0, which is more physically consistent with maximum erosion that takes 

482 place at that point. Based on the JET method, the critical shear stress is defined assuming 

483 a linear relationship between erosion rate and shear stress at the impingement point [13]. 

484 While Hanson and Cook [13] assumed maximum shear stress causes the maximum scour 

485 hole at the impingement point, the applied distribution of wall shear stress shows a zero 

486 value right at r/H = 0. Although the maximum shear stress based on the TKE distribution 
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487 is not at r/H = 0, it would partially justify the maximum erosion at that point by showing 

488 a high value dXr/H = 0 (Fig. 15). It would be useful to note that high mean pressure and 

489 also high fluctuations of pressure near the impingement plane especially at the impingement 

490 point could be considered as another driver of erosion [66]. Hence, the assumption of linear 

491 relationship between erosion rate and shear stress may not be valid at the impingement point. 

492 The position of the peak shear stress based on the Reynolds stress occurs at a greater 

493 distance from the jet centerline (impingement point) in comparison to the distribution based 

494 on TKE. Haehlen and Dade's [2] experiments determined the peak erosion rate location at 

495 about r/H = 0.1 for H/do = 154.2, similar to the location of peak wall shear stress at 

496 r/H = 0.14 reported by Beltaos and Rajartnam [22]. In this study, the peak shear stress 

497 occurred at r/H = 0.06 based on TKE and r/H = 0.30 based on Reynolds stress. The lower 

498 entrainment and lateral spreading rate of the free jet result in a thinner free jet in comparison 

499 to those formed in unconfined conditions. When a thinner jet impinges on the wall, the peak 

500 of the shear stress should occur at a smaller radial distance from the impingement point. This 

501 result agrees well with the wall shear stress based on TKE. There is good consistency between 

502 measured spreading rate of the free jet and position of the peak shear stress based on TKE. 

503 The spreading rate of the free jet is about 33 % less than previously reported values, and the 

504 wall shear stress peak, based on TKE, occurs about 33 % closer to the nozzle in comparison 

505 to the wall shear stress distribution of Phares et al. [23]. 

506 The peak values of the shear stress distributions found here are more than 2.4 times larger 

507 than those of Eq. (13). Values of C s were found to be 0.38 for both Runs 1 and 2, respectively, 

508 in comparison to the value of 0.16 proposed by Beltaos and Rajaratnam [22]. As previously 

509 discussed, their wall shear stress distribution appears theoretically invalid. The peaks of Run 

510 1 and Run 2 are, respectively, 28 and 53 % higher than the corresponding values in Phares 

511 et al. [23] distribution. It appears that jet confinement affects the shear stress distribution on 

512 the wall by decreasing the entrainment of fluid and by increasing the momentum transfer. 

513 This increased momentum accumulates in the relatively thin impinging jet and significantly 

514 increases the maximum shear stress exerted on the wall in comparison to the distributions 

515 predicted by Phares et al. [23] from the jet centerline up to a radial distance of about 0.4H 

516 (Fig. 16). The error bars on Fig. 16 are the 95 % confidence level for the extrapolation of TKE 

517 to the impingement wall. The confidence levels show that there is a meaningful difference 

518 between shear stress in the confined condition and all other relations for unconfined conditions 

519 for r/H < 0.4. 

520 There are several differences in the shear stress distribution based on TKE for Run 1 and 

521 Run 2 close to the impingement point. These differences might be due to the nozzle velocity. 

522 Since a wide range of nozzle velocities was not tested, no conclusion about the impact of 

523 this variable can be drawn. The distribution of shear stress in Run 2 is chosen for further 

524 analysis as its nozzle velocity was higher than for Run 1 and the distribution of shear stress 

525 appears better developed. To fit a curve through this distribution, the shear stress results based 

526 on TKE are divided into two zones: from r/H = 0 to 0.38, and r/H > 0.38 (Fig. 16). In 

527 the first zone, shear stress increases from a value of about 65 % of the maximum at the jet 

528 centerline to the maximum at r/H = 0.06, and it then decreases to a value of about 41 % of 

529 the maximum at r/H = 0.38. It is interesting to note that high values of shear stress between 

530 r/H = 0 and 0. 19 physically correspond to the lateral extent of the free jet since the free jet 

531 full- width right before entering the impingement region is slightly less than 5do (~ 0.20//; 

532 see Fig. 10) for the right side of the jet. A curve can be fitted through data between r/H = 0 

533 and 0.38 with R 2 = 0.99, 
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Fig. 16 TKE-based shear stress compared to literature shear stress distributions and best fit curve. The error 
bars shows the 95 % confidence level for the extrapolation of TKE to the impingment wall 
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(18) 



Between r/H = 0.19 and 0.38, the shear stress variation is low. In this zone, the radial 
velocity strongly decreases (see Fig. 7c) while the wall jet width is slowly increasing (see 
Fig. 7b), so it makes physical sense to have a low variation in the shear stress in this region. 
In the region r/H > 0.38, the shear stress distribution follows that of Poreh et al. [46], but 
with slightly higher values (about 13 % greater on average). 



540 4.3 Implications for sediment erodibility assessment 



541 As previously noted, the JET is an ASTM standard used for sediment erodibilty assessment in 

542 the field, and the results presented here address an important aspect of the analytical method 

543 developed by Hanson and Cook [13], which is the effect of confinement on bed shear stress 

544 at impingement. The method developed by Hanson and Cook is based on the presumption 

545 that the shear stress distribution of an impinging jet within the apparatus is identical to that 

546 observed in unconfined conditions. This study shows unequivocally that confinement can 

547 significantly alter the maximum value of this shear stress and its radial distribution from the 

548 impingement location. While the maximum shear stress for the conditions of this study and 

549 the particular confinement ratio was about 240 % higher than the value based on the JET 

550 methodology, the maximum shear stress is expected to be affected by the confinement ratio. 

551 Recently, a miniature version of the JET device ("mini" JET) was developed for easier use 

552 in field applications, and this device was compared to the original JET [29,67]. The "mini" 

553 JET is smaller and lighter than the original jet and has a smaller nozzle diameter (3.2 mm) 

554 and cylindrical tank diameter (about 0. 1 m) as compared to the original JET (6.4 mm and 0.3 

555 m, respectively). Note that the cylindrical tanks of the original and "mini" JET apparatuses 

556 are much smaller than the box used in the present study, so the effect of confinement is 

557 expected to be more severe for the JET. Al-Madhhachi et al. [29] compared the erodibility of 

558 carefully prepared sediment samples using both the original JET and the "mini" JET. They 

559 found that the soil's measured erodibility coefficient kd was the same between the original 
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and the "mini" JET, but the critical shear stress r c as determined using the "mini" JET was 
significantly lower (more than an order of magnitude lower) than that derived for the original 
JET. Al-Madhhachi et al. hypothesized that this difference in r c was due to the scales of the 
jet nozzles and the compaction method, and therefore introduced an adjustment coefficient to 
match the critical shear stress from the "mini" JET to the one from the original JET. Based on 
the results presented herein, it is more likely that the assumed shear stresses within both the 
original and "mini" JET devices on the sample surface are not equivalent. As noted above, 
the JET methodology under predicts the applied shear stress in the original apparatus due to 
the effects of confinement. The apparent discord between the original JET and the "mini" 
JET, i.e., the severe under prediction of the critical shear stress of the soil, likely is due to 
the stronger confinement conditions. Note that the confinement ratio of the "mini" JET is 
about 1,020, compared to the value of about 2,270 for the original JET. It is hypothesized 
here that the observed differences in the original and "mini" JET test results are due to the 
confinement effects. 

The size of the resultant scour hole created by the JET apparatus also might affect the shear 
stress distribution. Weidner [24] showed that for wide and shallow scour holes, the effect of 
scour hole shape on the shear stress distribution of the jet inside the hole is negligible. Thus, 
the maximum wall shear stress in these conditions can be found using the flat surface shear 
stress distribution (Eq. 18), 

(H e /d 0 r 

where C s is the shear stress coefficient and H e is the initial impingement height plus the scour 
hole depth. For the confined conditions presented here, C s is 0.38. The stronger confinement 
of the JET apparatuses is expected to increase wall shear stresses and increase values of C s 
in Eq. 19, with C s for the "mini" JET being greater than for the original JET. 



584 5 Conclusions 

585 The hydrodynamics of an axisymmetric impinging jet commonly employed to assess sedi- 

586 ment erodibility were investigated using PIV in a laboratory tank. The PIV apparatus was 

587 able to measure the velocity in both the near- and far-field regions of the free and wall jet, 

588 and these measurements agree well with the general behavior of impinging jets described 

589 in the literature. Mean and turbulent velocities follow self- similarity profiles in the free jet 

590 and wall jet region. Two major phenomena were observed in the confined conditions: (1) 

591 low entrainment because of limited source of fluid to be entrained into the jet, and (2) an 

592 increased rate of momentum transfer in the free jet because of the added momentum by the 

593 secondary flow. These two phenomena appear to play key roles in modifying the behavior of 

594 a confined jet in comparison to a jet in unconfined conditions. For example, the axial diffu- 

595 sion coefficient Cd was found to be a function of nozzle velocity and on average about 20 % 

596 higher than previously reported. The lateral spreading rate also was found to be a function 

597 of nozzle velocity and on average was 33 % lower than previously reported values. This 

598 low-spread jet with high diffusion in the axial direction altered the shear stress distribution 

599 at and downstream of the impingement point. 

eoo Wall shear stress was found based on the extrapolation of Reynolds stress and TKE. Using 

601 TKE to estimate shear stress showed a more consistent distribution as compared to other 

602 observations, particularly since it resulted in non-zero shear stress at the impingement point, 

603 which is typically where maximum erosion is observed. Based on TKE, the maximum wall 
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shear stress occurred 33 % closer to the impingement point, and with 240 % higher values, 
in comparison to the standard values recommended for use in the JET methodology. A semi- 
empirical equation (Eq. 19) was developed to predict the applied shear stress by the JET in 
the field. The value of shear stress coefficient C s in this equation should be reevaluated for the 
confined conditions of the original and "mini" JETs. Since these apparatuses are cylindrical 
and much smaller than the box in these experiments, the impact of jet confinement may be 
expected to be even more severe for the JETs. The present findings hypothesize confinement 
effects as the reason for observed inconsistencies between the "mini" and original JETs, and 
redefining the apparatus coefficients and relations should improve soil and sediment erosion 
assessment in the field. Further studies are needed to understand those effects in the actual 
conditions of the jet devices, to compute or measure the shear stress in the boundary layer 
with higher resolution to eliminate the errors induced by extrapolations, and to examine the 
role of pressure on erosion of sediment near the location of jet impingement. 
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